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Abstract: 

The relativistic corrections to the productions of double J/ip , double T, and J/ip + T at 
the Tevatron and the LHC were investigated within the frame of nonrelativistic QCD. The 
ratios of short distance coefficients between relativistic correction and leading order result for 
color singlet and color octet states in large pt limit are approximately —1,-11/3, respectively, 
for pair production. And for J/tp + T process the ratio is —11/6 for the CO channel. The K 
factors of relativistic corrections for color singlet and color octet double J/tp production are 
—0.23 and —0.84 with adopting v 2 = 0.23. The relativistic corrections significantly dilute the 
difference between the shape of color-singlet differential cross sections and the color-octet's at 
leading order. Also our results show that the single parton scattering model may be enough 
to explain the LHCb data for the J/tp pair production. 
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1 Introduction 

As the most achievable effective field theory to study production cross sections and decay 
rates for heavy quarkonium, the nonrelativistic QCD (NRQCD) [1] has been widely accepted. 
In the NRQCD approach, a heavy quarkonium can be described as a superposition of a 
short distance coefficients QQ n { 2s+1 Lj) mutiplied by a long-distance NRQCD matrix element 
< ( 2s+1 Lj) >, here, s, L, J denote the quantum numbers for the spin, the orbital angular 
momentum, the total angular momentum of the QQ intermediate state, respectively, and n 
stands for the color for the heavy-quark-antiquark pair state. The short distance coefficients 
indicate the creation or annihilation of a heavy quark pair can be calculated perturbatively 
with the expansions by the strong coupling constant a s . And the long-distance NRQCD 
matrix element is the probability of the intermediate heavy-quark-antiquark pair to evolve 
into an observed quarkonium nonperturbatively. In general, the intermediate state can be in 
a conventional color-singlet or a color-octet state. By introducing the color-octet mechanism, 
one can get the infrared-safe calculations for the decay rates of P-wave quarkonium states [2-9] 
and D-wave[10-12] decay widths. Besides, one may also explain the Tevatron data on the 
surplus production of J / tp and tp' at large pt by fitting the data and extracting the color-octet 
long-distance NRQCD matrix elements. 

Leaving aside the successes NRQCD have achieved, there are still puzzles in the com- 
parisons between the NRQCD predictions and the experimental data. The initial one is 
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the transversely polarized quarkonium prediction at the Tevatron failing in confronting with 
experimental measurements. This puzzle is called a "smoking gun" for the color-octet pro- 
duction mechanism. And recent ones are the large discrepancies between the leading order 
calculations of inclusive[13-15] and exclusive[16, 17] charmonia productions and the experi- 
mental data [18-20] at B factories. Now the most achievable efforts to resolve these puzzles 
seem to be higher order corrections. The next-to-leading order (NLO) QCD corrections [21- 
29], relativistic corrections [30-34], and 0(u s v 2 ) correction [35, 36] largely compensate the 
discrepancies between theoretical values and experimental measurements at B factory. The 
B factories analysis also show that the the color-octet matrix elements may be overestimated 
previously [24, 37, 38]. And the contributions of higher order corrections for charmonium 
prediction and polarization at hadronic colliders are still significant [37, 39-47]. There are 
more information in the Ref. [48, 49]. 

Most recently, the LHCb measured the J/ip pair production at a center-of-mass energy 
of y/s = 7TeV, with an integrated luminosity of 35.2 pb~ l and give the differential cross 
sections as a function of the J/ip pair invariant mass [50]. These predictions are thought to 
be useful to study the color-octet mechanism in quarkonium production by distinguishing the 
color-siglet contribution and color-octet contribution at different px region[51-55]. However, 
theoretical predictions within the single parton scattering mechanism are small. So double 
parton scattering mechanism was introduced in Ref. [56-58]. The predictions combining the 
double parton scattering and the single parton scattering lead a better agreement with data 
than that with the single parton scattering alone. 

However, as mentioned above, high order calculations including QCD and relativistic 
corrections to both CS and CO states are essential for heavy quarkonium production in 
the single parton scattering process. A complete next leading order calculations to double 
quarkonium production would be helpful to understand the multiple quarkonium production 
mechanism, COM mechanism, and the double parton scattering effects, at hadron colliders. 

Based on this motivation, in this paper, we studied the effects of relativistic corrections 
to both the CS channel and CO channel of production of double J/ip , production of double 
T at the LHC and the Tevatron. And we also discussed the relativistic corrections to the 
J/ip + T process for a comparison. 

This paper is organized as follows. In section II the NRQCD factorization and the 
standard methods to calculate the short distance coefficients of relativistic corrections are 
presented. In section III, we would give the numerical results and make a discussion. Finally 
there is a summary. 

2 The method of the calculations on relativistic corrections in NRQCD 

The differential cross section of two heavy quarkonium production can be discribed by con- 
voluting the cross parton level sections in the parton model as the following expression: 

da(p + p(p)^H 1 + H 2 + X) = j dx 1 dx 2 f a / p (x 1 )f b/p{p) (x2)da(a + b^H 1 + H 2 ), (2.1) 

a,b 
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where f a {b)/p(p){ x i) ls the parton distribution function(PDF), and Xj is the parton momentum 
fraction denoted the fraction parton carried from proton or antiproton. 

Within the NRQCD framework, the parton level cross sections can be divided into two 
parts: short distance coefficients and long distance matrix elements: [1] 

a{a + b^H, + H 2 ) = ^ f -^-4 ( Q l°m |0)(0|O^ |0>. (2.2) 

van Ql Q2 

On the right-hand side of the equation, the cross section is expanded to sensible fock states. 
F rnn noted by the subscript m, n, i.e. short distance coefficients, which can be obtained 
matching with the perturbative QCD, describe the process that produces intermediate QQ 
pairs in short range before harmonization to the physical meson state. (0|O^ n |0), the long 
distance matrix elements, represent the harmonization process that QQ evolutes to the CS 
final state by emitting soft gluons. The long matrix elements can be extracted from the 
experiment data[59, 60], determined by potential model[61] or lattice calculations [62-64]. 
0^ t n are local four fermion operators. The factor of nig" 1 '™ -4 is introduced to make F mn 
dimensionless. 

In this thesis, our calculations only refer to the S-wave spin-triplet quarkonium (J/ip,T) 
and its wave function looks like 

\H)=0(l)\QQ( 3 S?))+0(v)\QQ( 3 pf ] )g) + 0(v 2 )\QQ( 1 S^ 

(2.3) 

So the fock state expansion of Eq.(2.2) may involve different combinations of these states in 
above expressions. 

For the computations of short-distance coefficients, there are different processes con- 
tributing to the two meson production at parton level. In this paper, we just consider the 
gluon fusion process which is the dominate one of the parton level processes, i.e. g(ki) + 
9(^2)— >QiQi + Q2Q2^Hi(pi) + H 2 {p2)- Then the Lorentz-invariant Mandelstam variables 
at parton level are writen by (Here we omit the hats on the top of the s, t, u for simplification) 

s = (P1+P2? = (h + h) 2 , 
t = (h-Pi) 2 = (k 2 ~P2) 2 , 

u=(k 2 -p 1 ) 2 = {k 1 - P2 ) 2 , (2.4) 

and they satisfy 

s + t + u = M 2 Hl +M 2 H2 . (2.5) 

In the LO computation, the momenta of the heavy quark/antiquark in the intermediate 
QQ pair are set to half of the meson momentum. In an orthodox method to deal with the 
relativistic corrections, the amplitudes would be expanded by the rest relative momentum of 
the QQ around the quark/antiquark mass. In this paper, we also adhere to this strategy. We 
first define the momentum of heavy quark/antiquark as Pq/Pq- In a arbitrary frame, they 
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can be written as the following expressions [30, 65, 66]: 



PQ = -jP + ?> 

1 . . 

Pq = ^P~Q- (2-6) 

where p boosts from the momentum of the quarkoniu m in rest frame of the meson, etc. 
(2E q ,6) and q boosts from (0, q). Here, E q = ijm 2 + \q\ 2 is the rest energy of quark(antiquark). 

= itiqv is a lorentz scalar and v is the three relative velocity of quark or the antiquark 
in the rest frame. We can estimate v from the Gremm-Kapustin relation in weak-coupling 
regime[67, 68] 



v 2 = v\ = vl= ™ ocn W ■ (2-7) 



- 2m pole 



2m 



QCD 
Q 



where subscript "1" and "8" represents the pair of quark and antiquark is in a CS or CO 
state. And rriQ CD is the mass of quark in NRQCD actions and m P Q £ is the pole mass of 
quark. Taking Mqq, m P Q le , and rrig 00 as the input parameters, we can estimate the value 
of v. For example, in this paper, we select Mj/^ = 3.1GeV and m® CD = m pole = 1-39 , and 
then get v 2 ~ 0.23. And M T = 9AGeV and m® CD = m P b ° le = 4.57 ,and then get v 2 ~ 0.1. 
With the GK relation 

one may use the above values of v in estimation of the sizes of next leading order LDMEs. 

Taking advantage of the approximate relation Mh = 2E q and Eqs. (2.4,2.5), the Mandel- 
stam variables may be expanded by the relative momentum \q\. s is the beam energy and 
is \q\ independence. To expand t,u, we write they in the frame of gluon mass. Then the 
expressions expanding to next order may be given as 

tQq^, \q 2 \) = 2E 2 ql + 2E 2 q2 - S - - l\h ( s , (2E ql ) 2 , (2E q2 ) 2 ) 

= t(0, 0) + Ft{2m Q1 , 2m Q2 )\q 1 \ 2 + Ft{2m Q2 , 2m Q1 )\q 2 \ 2 + 0(v 4 ), 

uQq^, \q 2 \) = 2E 2 ql + 2E 2 q2 - | - l\i (s, (2E q2 ) 2 , (2E q2 ) 2 ) 

= u(0, 0) + Fu(2m Q1 , 2m Q2 )\q 1 \ 2 + Fu(2m Q2 ,2m Q1 )\q 2 \ 2 + 0(v*). (2.9) 

where A(x, y, z) = x 2 + y 2 + z 2 — 2{xy + xz + yz) and 

4{[t(0, 0) + u(0, 0)](t(0, 0) + y 2 ) - 2y 2 [t(0, 0) + x 2 }} 



Ft(x,y) 



s[s - 2{x 2 + y 2 )] + (x 2 - y 2 ) 2 ] 



_ 4{[t(0,0)+n(0,0)](n(0,0)+y 2 )-2^K0,0)+x 2 ]} 
tU ^' y )- s[s -2(x2 + yZ)} + (x2-y^} ' ^ iUj 
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For the case that the final two mesons are identical particals, the expanding expressions would 
be simplified to 

4[t(0,0) +4:171%] 

= t 0,0) - L 4- 2 qJ (|gi| 2 + |<f 2 | 2 ) + 0(« 4 , 

u(qi,q 2 ) = u(0, 0) - 1 2 0J (|gl| 2 + |g 2 | 2 ) + 0(v% (2.11) 

The short distance coefficients can be evaluated by matching the computations of per- 
turbative QCD and NRQCD: 

da =Y , F ™ n , A (O\O% 101 \O)(O\O% 2Q ~ 2 \O) . (2.12) 

pert QCD m Q\ m Q2 ' pert NRQCD V ' 

The covariant projection operator method should be adopted to compute the expression on the 
left-hand side of the equation. Using this method, spin-singlet and spin-triplet combinations 
of spinor bilinear in the amplitudes can be written in covariant form. For instance, we give 
the spin-triplet projection, 

E.,«000(S) & 'I h S M ) = ^ +mQ) (| - i - m Q )^^ + i + m Q ), (2.13) 

where e denotes the polarization vector of the spin-triplet state. Furthermore in our calcula- 
tion the Dirac spinors are normalized as uu = —vv = 2rriQ. Then the differential cross section 
satisfies: 

da(g + g^Q 1 Q l (^ +1 L [ jf) + Q 2 Q 2 ( 2s ^ 1 L [ j 2 8] )) 

+ 9^{QiQi){ 2si+1 L [ jf) + (Q 2 Q 2 )( 2s2+1 L [ ;f ] ))\ 2 
x(0|O Hl ( 2si+1 ^Ji)|0)(0|O H2 ( 2s2+1 L J2 )|0). (2.14) 
where the amplitude M is calculated by the following expression 
M(g + gMQiQi)( 2s+1 L [ }' 8] ) + (g 2 Q 2 )( 2s+1 4 1,81 )) 

= ^2 ^(o'* 1 ' 2 , ^ Jl,Jlz ^2 ,S2] 2'^ 2 ' J 2,^2 Z )(3«;3j|l,8a)(3/c;3l|l,86) 

S1S1S2S2 ijkl 

A(a + b^Q\ + Qi + Q\ + Q J 2 ). 

Next, we first expand the amplitude in powers of \q\\ and \q 2 \ to next order in v 2 (Here 
we only give the S-wave case and it is the only case using in our consequent calculations) 

M(g + g^ (QiQi) + (Q2Q2)) = (^) 1/2 (^) l/2 M( qi , q 2 , E ql , E q2 ) 

1 d 2 M 1 d 2 M 

= M(o,o,m Q1 ,m Q2 ) + (arfrf^j + ^rf^g 

d[(m Q1 /E ql )^M] d[(m Q2 /E q2 )^M] ^ n 
+ ^ 3E^ + ^ 8E^ XO.O.mgi.mga) 

+0(|gi| 4 ) + 0(\q 2 \ A ) = + -M (2) + 0(|gi| 4 ) + 0(|g 2 | 4 ). (2.15) 
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The factor (^) 1/2 comes from the relativistic normalization of the Fock state. For the 

Eg 

a P 

S-wave case, q a q^ = h\q\ 2 (— g a/3 + ^4r-) = i|g| 2 II Q ' 3 and the odd-power terms of q vanish. 
Subsequently, by multiplying the complex conjugate of the amplitude, the amplitude squared 
up to the next order can be obtained: 

\M\ 2 = M°M°* + (M°M 2 * + h.c.) + 0{v 4 ) 
= M°M°* 

+ 0(i; 4 ). (2.16) 

Form the expressions, one can see the relativistic corrections contain four parts to the two 
quarkonum production. To evaluate them, Eqs. (2.9 2.11) would be used. 

2.1 J/ip(T) pair production 

In Ref.[55, 69, 70], the production of J/ip pair was discussed in NRQCD. In the fock state 
expansion of the cross section, only two combinations give the dominate contributions and 
other ones are suppressed in all the pr region, i.e. psf^)! 3 ^ 1 ^) state and pS^g)! 3 ^ 8 ^). 
The former item gives dominate contribution at relative lower px region and the latter becomes 
important along with px increasing. In this paper, our computations restrict relativistic 
corrections of these two combinations. Therefore the differential cross section in Eq.(2.2) up 
to next order in v takes the following form 



da(g + g^J/ij + JM) = 1 1 , 1 ' (0\O J ^( 3 S^ ] )\0)(0\O J ^( 3 s\ l] )\0) 

+ 2 x G( bl / 1 j (O|P J ^( 3 M ] )l°)(O| 0j/V '( 3 M ] )|O) 



h 1 bl , bl J (O|0 J ^( 3 5[ 8] )|O)(O|0 J ^( 3 5i 8l )|O) 



^f3c[ 8 l3o[ 8 ]\ 

+ 2 x 1 \ 1 1 {0\V J/lp ( 3 sf ] )\0){0\O J/ ^( 3 sf ] )\0). (2.17) 
m c 

where F (...), G(...) are the short distance coefficients and for leading order and relativistic 
correction respectively. The factor 2 in the relativistic terms is given because the contributions 
of relativistic correction for identified particle are naturally equal. So we can calculate only 
the contribution of one J/tp, then multiply by 2 to get the total correction. The four-Fermion 
operators O h (A) with the dimension-6 and the dimension-8 four-Fermion operators V H (A) 
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for the relativistic correction operators are defined as[l] 



o H ( 3 s?) = V^x-xW 

H ( 3 sf ] ) = ^aT a X -X^T a i; 



(2.18) 



where % and i/j are the Pauli spinors describing anti quark creation and quark annihilation 
respectively, a is the Pauli matrices and D is the gauge-covariant derivative with ft = f>-£). 
T is the SU(3) color matrix. 

The above analysis is also adequate to the production of T pair. We use FeynArts[71] 
to generate Feynman diagrams and amplitudes and FeynCalc[72] handle the amplitudes. 
For double J/ip's and double T's process there are 31 Feynman diagrams for color-singlet 
channel and 72 Feynman diagrams for color-octet channel. We will give the results of the 
short distance coefficients to leading order and relativistic corrections in the Appendix. Our 
leading order results to color singlet and color octet channel are consistent with the ref [55]. 

2.2 J/ip + T channel 

As the discussion in Refs.[55, 70], the events to double J /ip under the luminosity of LHC is 70 
and it is difficult to study COM by use of the process of J/ip pair production for the LHC ex- 
periments as well as for the di-T pair production. They study the process pp — > J/tp + T + X 
and find it is nice probe to COM with enough events at LHC. To study relativistic cor- 
rections to this process is a meaningful work. At leading order of 0(a s ), the CS channel 
cc( 3 sj 1] ) + bb( 3 SP) have no contributions. The CO channels include cc( 3 S[ 1] ) + bb( 3 Sf ] ), 
cc( 3 sf ] ) + bb( 3 sP) and cc( 3 sf ] ) + bb( 3 sf ] ) combinations. In our calculations, we neglect the 
relativistic corrections to T production which is negligible compared with the corrections to 
J/ip. Therefore we have the obvious expression to next order in v 



da(g + g -> J/if> + T) 



+ 



+ 



+ 



+ 



+ 




(0|O J/,/, ( 3 5f ] )|0)(0|C' T ( 3 5j 11 )|0) 



(0|7 ?J/ ^( 3 5j 8] )|0)(0|O T ( 3 5| 11 )|0> 



(O|0 J/ ^( 3 sf ] )|O>(O|0 T ( 3 sf ] )|O) 



(0|P J/,/ '( 3 5j 11 )|0)(0|(!7 T ( 3 sf 1 )|0) 



(0\V J/ ^( 3 sf ] )\0){0\O r ( 3 sf ] )\0). (2.19) 



(0|O^( 3 S , [ 1] )|0)<0|C x ( 3 S'l 8] )|0) 
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The operators are defined in Eq.(2.18) and the results of the short distance coefficients can 
also be found in the Appendix. 



3 The Numerical Result and Discussion 



Using the method of phase space integration in Ref. [59], we can carry out the numerical 
calculation the cross section at the LHC and Tevatron respectively. Here we ignore the 
relativistic correction come from the phase space. Fortran program was used to calculate 
the integration. We choose input parameters: the mass of charm quark m c = l.bGeV and 
bottom quark m& = A.lGeV; the factorization scale \i taken as transverse mass \x = itlt = 
(ArriQ 2 +pt 2 ) j specifically we use the parton distribution function CTEQ6L1[73]. The running 
coupling const a s is evaluated by the formula of CTEQ6L1. We employ the following values 
for the NRQCD matrix elements for J/V> and T [61, 74-76]: 

(0|O J /^( 3 5| 1] )|0) = l.3GeV 3 , 

(0\O J/tp ( 3 sf ] )\0) = 3.9xl(T 3 GeU 3 . (3.1) 



(O|0 T ( 3 S[ 1] )|O) 
(0|O T ( 3 5f 



9.2GeV 3 , 
l.bxW^GeV 3 . 

As discussion in Sec. 2, we use Eq.(2.8) to estimate the relativistic LDMES. 



(3.2) 



3.1 J/ip(T) pair production 

Under the large pt limit approximation, 



2 2 2 

m m m 

— ~ < — k- ~ 

u t p T 



(3.3) 



where m is the mass of J /tp or T, we can expand the short distance coefficients with m. The 
ratios of the short distance coefficient between LO (F) and its relativistic correction (G) in 
lowest order in the expansion are: 



R( 3 S? 3 S?) 



R( 3 sf ]3 S?) 



G( 3 S[ LiS S 



p T >m 



F( 3 S[ 1]3 S[ 



G( 3 sf ]3 sf ] ) 



F (3 S fh S fh 



1 

11 
Y 



(3.4) 



The ratio of 3 sf^ is approximately —11/3 this ratio of color octet state is consistent with 
Ref. [68, 77]. In large px limit, the dominate contribution of this subprocess is g* — > cc 
The propagator of virtual gluon g* is proportional to 1/E 2 This term offers a factor of —2 
to the ratio i?( 3 s[ 8] ). 



- 8 - 



In Fig.l, we present the ratios at the Tevatron with y/s = 1.96TeV and at the LHC with 
y/s = 7TeV or y/s = lATeV of J/ip and the curves are very close at large px and consistent 
with above analysis. In small px region, for CS, corrections are less than the values in large px 
limit. On contrary, there is a sharp suppress in small px region for CO. That can obviously be 
shown in Fig. 2, etc. differential cross sections of px- The LO cross sections to CS and CO state 
suppressed by factor of 23% and 84% when involving the effects of relativistic corrections. 
For the LO case, the CS and CO curves cross at about px = 14GeV . As px increasing the 
contributions to CO state would surpass that to CS. And this provide a probe to check the 
COM experimentally. But in our calculations as shown in the diagrams, relativistic corrections 
reduce the CO contributions largely and make the difference not apparently between the CSM 
results and the COM results. Furthermore, we make a prediction of the cross sections at the 
LHCb with y/s = 8TeV as shown in Fig. 3, COM with relativistic corrections contribute little 
to total cross sections that is similar with the cases mentioned above. For total cross sections, 
COM effects are also not obviously, and this can be seen in Tab.l. 
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Figure 1. The ratios of the short distance coefficient between LO F and its relativistic correction G 
of p + p — > J/tp + J/ip process at the Tevatron with y/s = 1.96TeV and at the LHC with y/s = 7TeV 
of = UTeV. 
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Figure 3. The predicted p + p — >• J/ip+ J/ip process pt distribution of deferential cross section at 
the LHC with ^/s = 8TeV of \y\ < 2.4 and 2.5 < y < 4.5 , respectively. 



In comparison of process of single J/tjj production in which QCD corrections are essential 
than relativistic corrections to CO states and give large contributions to total cross sections. 
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3 SiW LO 


3 SiW RC 


3 S^1 LO 


RC 


total cross section 


(GeV) 


(nb) 


(Pb) 


(Pb) 


(Pb) 


(nb) 


1.60 


5.18 


-13.90 


1.63 


-1.94 


5.16 


1.55 


6.92 


-22.12 


2.20 


-2.61 


6.89 


1.50 


9.31 


-34.68 


3.00 


-3.55 


9.28 


1.45 


12.65 


-53.09 


4.13 


-4.88 


12.60 


1.40 


17.34 


-80.96 


5.73 


-6.77 


17.26 



Table 1. The total cross sections for color-singlet channel, color-octet channel and their relativistic 
corrections at LHC with y/s = 7TeV : integrating over 2 < y < 4.5 and pr < WGeV and we choose 
charm quark mass for input parameter from lAGeV to l.QGeV. "RC" at the head of the table denotes 
the relativistic contributions. 



It is expected for QCD corrections to CS and CO states to give a large contributions. 

To investigate the effects of relativistic corrections comparing with experiment, we give 
the differential cross sections as functions of the invariant mass of double J/tp and of the 
rapidity in Fig. 4 to pair production with the same constraints as the LHCb, i.e. 2 < y < 4.5 
and ^/s = 7TeV. From the figures, one can see the main part of the theoretical calculations 
concentrates in the low region of the invariant mass with a sharp peak at about 6.5GeV , but 
LHCb data shows a peak at 7 ~ 8GeV. And theoretical value of the peak is nearby 5.5nb/GeV 
which is two times large than the experimental peak value. In the low region less than about 
7.5GeV, relativistic effects show a positive correction to LO curve with a maximal enhanced 
factor 1.08. In the region large than 7.5GeV, relativistic effects suppress the LO curve but 
not obviously. Here we only consider the contributions to direct pair production and the 
feed-down process from excited state will contribute nearly 30% to direct production rate [56]. 
Taking the feed-down corrections into account, the theoretical prediction to differential cross 
section within the invariant mass 6GeV ~ 8.5GeV would be large than the experimental 
values. This situation is very different from the calculations in Refs.[56, 57]. They showed 
that the single parton scattering contributions were too small to fit the experimental data and 
introduced double parton scattering mechanism. In our calculations, only LO contributions 
would be large enough comparing with the experimental data. But as they discussed, single 
parton scattering result peaks at too low a invariant mass value and with a too large value, and 
the relativistic corrections do not relieve the puzzle but enhances the peak value. Recently, 
an soft-collinear effective theory(SCET) has been developed and seem to provide a possibility 
to suppress the end point peak considerably and move it from the end point[78]. Besides, 
may NLO QCD redia corrections give a large negative contribution especially for the peak? 
It is expected the problem would be solved when involving these corrections. 
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Figure 4. Left figure is the differential cross section to double J/ip production as a function of the 
invariant mass tuj/^j/^ at LHC with y/s = 7TeV and 2 < y < 4.5, and dash and solid lines correspond 
to LO and NLO(u 2 ) results, respectively. Right figure is for rapidity distribution, dash and solid lines 
also correspond to LO and NLO(i; 2 ), respectively. In both two figures, it were set to Mj/^ = 3.1. 




Figure 5. The figure is about the differential cross section to double J/ip production over rapidities 
of J/ip J/ip at LHC with y/s = 7TeV and 2 < y < 4. 5, and yfs = 8TeV and 2.5 < y < 4.5 respectively. 



3.2 J/ip + T channel 



We give the differential cross sections as a function of pt at the LHC with \rj\ < 2.4 when 
choosing the beam energy of the LHC and the LHCb, \fs = 7TeV and \fs = 14TeV, respec- 
tively. 

As discussed in Sec. 2 we consider only the relativistic corrections of the J ftp part. The 
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p T (GeV) p T (GeV) 



Figure 6. The deferential cross section ofp + p— > T + T process as a function of pt at the LHC with 
\r)\ < 2.4 and y/s = 7TeV , y/s = UTeV, respectively. 



radios of F/G in the large px limit are given below 



G 


cc( 3 S? 


,bb( 3 sf ] )\ 




11 


F 




,bb( 3 S [ 8 ] )\ 






G 




M 3 s?)] 




11 


F 




,bb( 3 S?)} 


p T >m 


~~6 



R[cc( 3 sP),bb( 3 S?)} 



R[cc( 3 sf ] ),bb( 3 S?)] 

(3.5) 

For the channel cc( 3 S^) + 66( 3 5[ 8 ^), the ratio is not a const result and its value approximate 
1.5 at the large pt region. In Fig. 7, we give the Pt distribution of the LO cross sections 
and the next leading order ones for different channels. The relativistic corrections for the 
production of J/tp + T are limited to change the LO CO contributions. This process may be 
still absorbing to research to be a probe to COM. 



4 Summary 

We calculate the relativistic corrections to the process p + p(p) — > J/ip + J/tjj + X ,p + p(p) — > 
T + T + X, and p + pip) ->• J/ip + T + X at the Tevatron and at the LHC based on NRQCD. 
The ratios of short distance coefficients between relativistic corrections and LO results for CS 
and CO states are approximately — 1, — in large pt limit for the J/ip and T pair produc- 
tions. And for J/ip + T the ratio is — -g- for the CO channel. 

If fix the LO LDMEs and estimated the NLO(v 2 ) LDMEs through the velocity scaling 
rule with adopting v 2 = 0.23, the differential cross sections are suppressed by 23% and 
84% to LO cross sections for CS channel and CO channel at large pt region, respectively. 
The relativistic corrections significantly dilute the difference between the shape of color- 
singlet differential cross sections and the color-octet's at leading order. So, it is difficult to 
study the color-octet mechanism in quarkonium production by distinguishing the color-siglet 
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Figure 7. The differential cross section of p + p — > J/ip + T process as a function of pt at the LHC 
with Vs = 7TeV and = UTeV with |r?| < 2.4 . 

contribution and color-octet contribution at different pt regime, if the relativistic corrections 
are considered in pair heavy quarkonia productions at hadron collider. 

We also compute the differential cross sections as a function of the double J/tp invariant 
mass and compare it with the recent LHCb experiments. Our results show that the single 
parton scattering contributions are large enough to fit experimental data and it would leave 
limit room to double parton scattering mechanism. 

Note: when this paper is preparing, another two thesises appeared which discussed the 
relativistic corrections of the pair production of J/tp and 7] c at the hadron collider [79, 80] but 
they were computed within the relativistic quark model. Comparing with their results, the 
corrections under the NRQCD frame are much smaller than that given by relativistic quark 
model. 
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5 Appendix 

5.1 short distance coefficients of g + g^-J/ip + J/ip and g + g — > T + T 

The LO short distance coefficients and their relativistic corrections ones for parton-level 
process g{k\) + g(k 2 ) ->■ J/ip(pi) + Jj^ip-i) and g(ki) + g(k 2 ) — >■ T(pi) + "f (p 2 ) are given in 
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this appendix. Here m = 2m c or 2m& represents the J/ip mass or T mass. (Here we adopt 
t(u) to represent i(0)(«(0)) directly for simplification) 



The LO short distance coefficient for CS channel is: 



F( 3 <?1 1]3 S[ 1] ) = „ , s7T \a/ 5 [2680m 24 - 14984m 22 (i + u) 

v 1 17 6561m 2 (m 2 -t) 4 (m 2 - u) A (-2m 2 + 1 + u) 8 I v ; 

+m 20 (31406t 2 + 89948tn + 31406u 2 ) - 16m 18 (1989t 3 + 12661t 2 u 

+1266Uu 2 + 1989n 3 ) + 4m 16 (4417t 4 + 57140t 3 n + 117714iV + 57140tn 3 + 4417-u 4 ) 

-4m 14 (1793i 5 + 38340t 4 u + 134119iV + 134119i 2 u 3 + 38340tu 4 + 1793u 5 ) + m 12 (2956£ 6 

+76406t 5 w + 361624tV + 571900tV + 361624tV + 76406tti 5 + 2956u 6 ) - 2m 10 (397t 7 

+1539U 6 u + 91227iV + 167593i 4 tt 3 + 167593tV + 91227iV + 15391tu 6 + 397u 7 ) 

+m 8 (47t 8 + 7642i 7 n + 73146t 6 u 2 + 150334t 5 u 3 + 132502i 4 u 4 + 150334tV + 73146t 2 n 6 

+7642tn 7 + 47u 8 ) + 2m 6 (10t 9 - 411t 8 n - 8951iV - 29063t 6 n 3 - 17653i 5 u 4 - 17653t 4 n 5 

-29063i 3 u 6 - 8951fV - 411tn 8 + 10u 9 ) + m 4 (t 10 - 66t 9 u + 2469tV + 12874iV 

+11928iV + 1164tV + 11928iV + 12874iV + 2469iV - 66tn 9 + u 10 ) + Am 2 t 2 u 2 (9t' 

-586t 6 u - 37tV - 394iV - 394tV - 37tV - 586iu 6 + 9u 7 ) + 2t 4 n 4 (349t 4 - 908t 3 u 

+1374iV - 908tu 3 + 349u 4 )l (5.1 
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The LO short distance coefficient for CO channel is: 



~ 972m 6 (m 2 - t) 4 (m 2 - uf (-2m 2 + t + uf 
- 104624m 28 + 522544?n 26 (i + u) - 2m 24 (1377 Alt 2 + 1836350tu + 137741u 2 ) 

-16m 22 (80955t 3 - 345827t 2 n - 345827i« 2 + 80955u 3 ) + 4m 20 t(667342t 4 - 701017t J 
-1566126tV - 701017tn 3 + 667342u 4 ) - 8m 18 (316747t 5 + 1677t 4 u + 195158t 3 u 2 
+195158iV + 1677tti 4 + 316747u 5 ) + m 16 (1455637t 6 + 281021i 5 u + 5910067t 4 n 2 
+14296294t 3 n 3 + 5910067tV + 281021tu 5 + 1455637u 6 ) - 2m 14 (270334t 7 - 179453t 6 n 
+2260581tV + 8407066tV + 8407066tV + 226058UV - 179453^ 6 + 270334?/) 
+2m 12 (66670t 8 - 239050t 7 ti + 1059585iV + 5792306iV + 7997354tV + 5792306t 3 u 5 
+1059585t 2 n 6 - 239050tu 7 + 66670u 8 ) - 2m 10 (10661i 9 - 110448t 8 tt + 408662tV 
+3075158t 6 n 3 + 5105335t 5 it 4 + 5105335tV + 3075158i 3 n 6 + 408662iV - 110448^ 8 
+10661ti 9 ) + 2m 8 (935t 10 - 25638t 9 n + 178326tV + 1288892tV + 2655522iV 
+3136494tV + 2655522t 4 it 6 + 1288892iV + 178326t 2 u 8 - 25638tii 9 + 935n 10 ) 
+2m 6 tu(2616t 9 - 65577t 8 n - 431588tV - 1010108£ 6 ii 3 - 1490351t 5 n 4 - 1490351t 4 ti 5 
-1010108t 3 n 6 - 431588tV - 65577iu 8 + 2616ti 9 ) + m 4 to(-243i 10 + 33354t 9 n 
+192771iV + 555338tV + 879860i 6 u 4 + 1154016i 5 ti 5 + 879860t 4 ti 6 + 555338t 3 n 7 
+19277UV + 33354tn 9 - 243ii 10 ) - 162m 2 t 2 u 2 (t + n) 3 (27i 6 + 88i 5 u + 174t 4 u 2 
+208tV + 174iV + 88tu 5 + 27u 6 ) + 243fcV(t + uf{t 2 + tu + u 2 f] (5.2) 
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The short distance coefficient to relativistic correction for CS channel is: 



{ 1 l> 19683m 4 (m 2 -t) 5 (m 2 --u) 5 (-2m 2 + t + u) 10 (2m 2 + t + -u) 
370752m 34 - 2181536(t + u)m 32 + 64[73726(t 2 + u 2 ) + 211649ut]m 30 

-16[244805(t 3 + u 3 ) + 1965309ut(t + u)]m 28 - 4[316463(t 4 + u A ) 
-8160588ni(t 2 + u 2 ) - 18782102u 2 i 2 ] ?n 26 + [5600590(t 5 + u 5 ) 
-9749802ui(i 3 + u 3 ) - 79767876u 2 i 2 (t + u)]m 24 - 2[2766605(t 6 + u 6 ) 
+5407076ui(i 4 + u A ) - 15844965u 2 £ 2 (t 2 + u 2 ) - 35085112u 3 i 3 ] m 22 
+2[l563020(t 7 + u 7 ) + 6369573nt(t 5 + u 5 ) + 2696929u 2 i 2 (i 3 + u 3 ) 
-456546u 3 t 3 (t + u)]m 20 - 4[295064(t 8 + u 8 ) + 1552733tit(t 6 + u 6 ) 
+1230476n 2 t 2 (t 4 + u A ) + 7465803u 3 t 3 (t 2 + u 2 ) + 17306424u 4 t 4 ] m 18 
+2[l57467(t 9 + u 9 ) + 858538-ui(i 7 + u 7 ) - 1932037 u 2 t 2 (t 5 + u 5 ) 
+3904293u 3 t 3 (t 3 + u 3 ) + 38310315ti 4 t 4 (t + u)]m 16 - 2[30250(t 10 + u 10 ) 
+134737-ui(i 8 + u 8 ) - 2493940n 2 t 2 (i 6 + u 6 ) - 5184812u 3 t 3 (t 4 + u 4 ) 
+15398042<u / 4 4 (t 2 + u 2 ) + 35880758tt 5 £ 5 ] m 14 + [8417(t n + u 11 ) 
+34279ui(i 9 + u 9 ) - 2497569u 2 t 2 (t 7 + u 7 ) - 10315939n 3 t 3 (t 5 + u 5 ) 
+1057162n 4 t 4 (t 3 + u 3 ) + 31198738u 5 t 5 (t + u)]m 12 + [ - 817(i 12 + u 12 ) 
-13616nt(t 10 + u 10 ) + 672840n 2 i 2 (t 8 + u 8 ) + 4910632n 3 t 3 (t 6 + u 6 ) 
+5008157n 4 t 4 (t 4 + u 4 ) - 7403400u 5 t 5 (£ 2 + u 2 ) - 13064072n 6 t 6 ] m w 
- [2l(t 13 + u 13 ) - 4794ut(t n + u 11 ) + 84325n 2 t 2 (t 9 + u 9 ) 
+1346073n 3 t 3 (t 7 + u 7 ) + 3171602n 4 t 4 (i 5 + u 5 ) - 481617n 5 t 5 (t 3 + u 3 ) 
-3414586n 6 t 6 (t + u)] m 8 + [7(i 14 + u lA ) - 2A8ut{t 12 + u 12 ) 
-2777u 2 t 2 (t w + u 10 ) + 197002u 3 t 3 (t 8 + u 8 ) + 973759u 4 t 4 (t 6 + u 6 ) 
+633832w 5 t 5 (i 4 + u 4 ) - 734333n 6 t 6 (t 2 + u 2 ) - 747796u 7 t 7 ]m 6 
+tu[{t 13 + u 13 ) + 137ut(t u + u 11 ) - 1320u 2 t 2 (i 9 + u 9 ) 
-171964u 3 t 3 (t 7 + u 7 ) - 249123n 4 t 4 (t 5 + u 5 ) + 36693n 5 t 5 (t 3 + u 3 ) 
+125000U 6 t 6 (t + u)] m A + 2t 3 u 3 [2(i 10 + u 10 ) + 761ut(t 8 + u 8 ) 
+36064u 2 i 2 (i 6 + u 6 ) - 14592u 3 t 3 (t 4 + u A ) + 9630u 4 t 4 (t 2 + u 2 ) 
-5362u 5 t 5 ] m 2 - 2iV [227(i 7 + u 7 ) + 6021ut(t 5 + u 5 ) 
-9217u 2 t 2 (t 3 + u 3 ) + 6041u 3 t 3 (t + u)] 



The short distance coefficient to relativistic correction for CO channel is: 

4^3 



a s ir 



m 34 



/-<C3o[8]3c[8K _ 

1 1 ; 729m 8 s 2 (m 2 - t) 5 (m 2 - uf{-2m 2 + t + u) 8 (2m 2 + t + u) 
- 7840128m 38 + 580182400 + u)m 36 - 16 [4351865(i 2 + u 2 ) + 36821254ut 
-8[22135741(t 3 + u 3 ) - 183792021tit(t + u)]m 32 + 4[l73739553(t 4 + u 4 ) 
-435034596ui(i 2 + u 2 ) - 916724794u 2 t 2 ] m 30 - 2 [576450337(t 5 + u 5 ) 
-575408499ui(i 3 + u 3 ) - 2090528238u 2 t 2 (t + u)]m 28 + 2 [612519223(t 6 + u 6 ) 
-216110060nt(t 4 + u 4 ) - 1597249863ii 2 t 2 (t 2 + u 2 ) - 1708173704n 3 t 3 ] m 26 
+2[ - 462600109(t 7 + u 7 ) + 86546766ut(t 5 + u 5 ) + 1187377312n 2 t 2 (t 3 + u 3 ) 
+911415039u 3 t 3 (t + u)]m 24 + [519274693(t 8 + u 8 ) - 285877238nt(t 6 + u 6 ) 
-2090011040u 2 i 2 (i 4 + u 4 ) - 1710759498u 3 t 3 (t 2 + u 2 ) - 1356254154n 4 t 4 ] m 22 
+ [- 221652624(t 9 + u 9 ) + 375759263nt(t 7 + u 7 ) + 1820975245u 2 t 2 (t 5 + u 5 ) 
+1804776207u 3 i 3 (i 3 + u 3 ) + 2287765845u 4 t 4 (t + u)]m 20 + 2[36158642(i 10 + u w ) 
-146055034«t(t 8 + u 8 ) - 613461077n 2 t 2 (t 6 + u 6 ) - 576102238u 3 t 3 (t 4 + u 4 ) 
-723890837n 4 t 4 (t 2 + u 2 ) - 1150404320n 5 t 5 ] m 18 + [ - 17656681(t n + u 11 ) 
+148210450ni(t 9 + u 9 ) + 547766115n 2 t 2 (t 7 + u 7 ) + 168319214ti 3 t 3 (t 5 + u 5 ) 
-601363505u 4 t 4 (t 3 + u 3 ) - 515358665u 5 i 5 (i + u)]m w + 2 [l515472(i 12 + u 12 ) 
-25773835ut(t 10 + u 10 ) - 63896082n 2 t 2 (i 8 + u 8 ) + 200994011u 3 t 3 (t 6 + u 6 ) 
+902427166n 4 t 4 (t 4 + u 4 ) + 1542963240n 5 t 5 (t 2 + u 2 ) + 1808885720u 6 t 6 ] m 14 
-[320244(t 13 + u 13 ) - 12358347ui(i n + u 11 ) + 10043059u 2 t 2 (t 9 + u 9 ) 
+412260621n 3 t 3 (t 7 + u 7 ) + 1621412381n 4 t 4 (t 5 + u 5 ) + 3279793470m 5 t 5 (i 3 + u 3 ) 
+4552155452u 6 t 6 (i + u)]m 12 + [l4662(t 14 + u 14 ) - 1965047nt(t 12 + u 12 ) 
+19157992-u 2 t 2 (t 10 + u 10 ) + 213718612n 3 t 3 (t 8 + u 8 ) + 861766390w 4 t 4 (t 6 + u 6 ) 
+1999127983u 5 t 5 (t 4 + u 4 ) + 3170423452m 6 t 6 (t 2 + u 2 ) + 3695142248u 7 i 7 ] m 10 
+tn[l90300(t 13 + u 13 ) - 7050997ni(i n + u 11 ) - 69336231u 2 t 2 (t 9 + u 9 ) 
-300672607n 3 t 3 (t 7 + u 7 ) - 7886428 llw 4 ^ 5 + u 5 ) - 1422124344-u 5 i 5 (i 3 + u 3 ) 
- 1887768238m 6 i 6 (i + u)]m 8 + tu[- 8019(i 14 + u 14 ) + 1359759ut(t 12 + u 12 ) 
+14495197ti 2 t 2 (t 10 + u 10 ) + 69109969w 3 t 3 (t 8 + u 8 ) + 205737389n 4 t 4 (t 6 + u 6 ) 
+417998865n 5 t 5 (t 4 + u 4 ) + 630439593u 6 i 6 (i 2 + u 2 ) + 718419166n 7 t 7 ] m 6 
-t 2 u 2 [I43856(t 13 + u 13 ) + 1824408iit(t 11 + u 11 ) + 10108398n 2 t 2 (t 9 + u 9 ) 
+33845791u 3 t 3 (t 7 + u 7 ) + 78611535n 4 t 4 (t 5 + u 5 ) + 1 32750949m 5 i 5 (> 3 + u 3 ) 
+172324471u 6 t 6 (t + u)) m 4 + 81t 2 u 2 (t + uf [81(t 8 + u 8 ) + 999*ut(t 6 + n 6 ) 
+2910u 2 t 2 (t 4 + u 4 ) + 5459u 3 t 3 (t 2 + u 2 ) + 6452u 4 t 4 ]m 2 

-2673t 3 u 3 (t + u) 7 {t 2 + ut + u 2 ) 3 ] (5.4) 
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5.2 short distance coefficients ofg + g^-J/ip + T 

The LO short distance coefficient for g + g — > cc( 3 <sj 8 ') + bb( 3 S^) channel is: 

F(cc( 3 S [8] )bb( 3 S [1] )) = 10a ^ 3 

{ { 1 ' { 1 )] 2A3m 3 M{M 2 - t) 2 (M 2 - u) 2 {-2m 2 + i + u) 2 s 2 

2m 8 M 2 + 2m 6 (7M 4 - 5M 2 (t + u) + tu) + 2m 4 (7M 6 - l3M 4 (t + u) 

+M 2 (7t 2 + 10tu + 7-u 2 ) - tit(t + «)) + m 2 (M 8 - 8M 6 (t + u) 
+M 4 (15t 2 + 16tu + 15u 2 ) - 4M 2 (2i 3 + 3t 2 u + 3iu 2 + 2u 3 ) + t 2 u 2 ) 
+2M 2 (M 4 (t 2 + tu + u 2 ) - M 2 (2t 3 + 3t 2 u + 3tu 2 + 2u 3 ) 

+ {t 2 + tu + u 2 ) 2 )}, (5.5) 

where m = mj/^, and M = my. 

The relativistic of short distance coefficient for g + g — > cc( 3 5[ 8 ') + 66( 3 S'[ 1 ') channel is: 

G(cc( 3 s[ 8] )&6( 3 S[ 1] )) = -20a 4 vr 3 /(729m 5 Ms 2 (M 2 - t) 3 (M 2 - u) 3 (4m 2 M 2 - (t + u) 2 )) 

/(-2m 2 + t + u) 3 x 96M 4 (2M 2 - t - u)m 14 + 16M 2 [61M 6 - 82(t + u)M 4 + (30(t 2 + u 2 ) 

+37ut)M 2 + 3tu(t + u)]m 12 + 4[52M 10 - 374(t + u)M 8 + (549(t 2 + u 2 ) + 290ut)M 6 
+(-227(t 3 + u 3 ) + 61ut(u + i))M 4 - tu(85(t 2 + u 2 ) + 2Uut)M 2 + 12t 2 u 2 (t + u)] m 10 
-2[424M 12 - 916(4 + u)Af 10 + 2(77(4 2 + u 2 ) + 934trf)M 8 + (867(i 3 + u 3 ) 
-U39ut(t + u))M 6 + (-505(4 4 + u 4 ) + 728ut(t 2 + u 2 ) + \778u 2 t 2 )M 4 - tu(241(t 3 + u 3 ) 
+735ut(u + t))M 2 + 2Qt 2 u 2 (t + u) 2 ]m 8 - 2[44M 14 - 512(4 + u)M 12 + 2(545(t 2 + u 2 ) 
+972ut)M 10 - (451(£ 3 + u 3 ) + 2589ni(i + u))M 8 + (-572(4 4 + u 4 ) + 1158ui(i 2 + u 2 ) 
+2264<u 2 i 2 )M 6 + (401(4 5 + ti 5 ) - 155ui(i 3 + u 3 ) - Wl8u 2 t 2 (t + u))M 4 + 2iu(41(4 4 + u 4 ) 
+128nt(t 2 + u 2 ) + 196u 2 4 2 )M 2 - 3t 2 u 2 (t + uj 3 ]m 6 + 2[l0(t + u)M 14 - (239(t 2 + u 2 ) 
+390ut)M 12 + 2(252(4 3 + u 3 ) + 611«t(* + u))M 10 - (117(t 4 + u 4 ) + 1676u 2 4 2 
+1027ut(u 2 + t 2 ))M 8 - (383(4 5 + u 5 ) + 2Uut(t 3 + u 3 ) - I57u 2 t 2 (t + u))M 6 
+(225(i 6 + u 6 ) + 485ut(t 4 + u 4 ) + 670u 2 4 2 (i 2 + u 2 ) + 684u 3 t 3 )M 4 - tu(56(t 5 + u 5 ) 
+205nt(t 3 + u 3 ) + 325u 2 t 2 (t + u))M 2 - t 2 u 2 (t + u) 2 (t 2 - 9ut + u 2 )]m 4 + [(71(t 3 + u 3 ) 
+I73ut{t + u))M 12 - (t + u) 2 (79(i 2 + u 2 ) + 90ut)M 10 - (211(t 5 + u 5 ) + 590ut(t 3 + u 3 ) 
+859u 2 t 2 (t + u))M 8 + {t + uf{37h{t 4 + u 4 ) + 634tit(t 2 + u 2 ) + 1006n 2 t 2 )Af 6 
-(156(£ 7 + u 7 ) + 827ut(t 5 + u 5 ) + 2082u 2 t 2 (t 3 + u 3 ) + 3115u 3 t 3 (t + u))M 4 
+tu(t + -u) 2 (108(i 4 + u 4 ) + 241ut(t 2 + u 2 ) + 33Au 2 t 2 )M 2 - 5t 3 u 3 (t + u) 3 ]m 2 
-22M 2 {M 2 - t){M 2 - u)(t + u) 3 [(t 2 + ut + u 2 )M 4 - (2(i 3 + u 3 ) + 3ut(t + u))M 2 



+(t 2 + ut + u 2 ) 21 



(5.6) 
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The short distance coefficient of gg — > cc( 3 S^) + bb( 3 S^) can be get through Exchanging 
M and m. The relativistic correction of this cross section is : 

^3 Q [W3o[8]^ 20a 4 7r 3 /(729m 3 M 3 S 2 ) 



G(cc( 3 S[ Ll )bb( 3 SP)) 



t) 3 (m 2 - u) 3 (-2M 2 + t + u) 4 {Am 2 M 2 - (t + u) 2 ) 



(16M 4 [14M 2 - 9(t + u)] m 16 + 4M 2 [300M 6 - 452(t + <u)M 4 + (161(t 2 + u 2 ) + 274ui)M 2 

-9(i 3 + u 3 ) + 13tu(< + u)]m 14 - [288A/ 10 + 1648(t + u)M 8 - 4(575(i 2 + u 2 ) + 894nt)M 6 
+16(32(i 3 + u 3 ) + 3ut{t + u))M 4 + (t + u) 2 (3(t 2 + u 2 ) + 574ui)M 2 - 4(< + ti) 3 (3(i 2 + u 2 ) 
+2ut)] ?n 12 - [800M 12 - 2080(t + u)M 10 + 16(39(i 2 + u 2 ) + 73ut)M 8 + 4(121(i 3 + u 3 ) 
-65ut(t + u))M e + (34(t 4 + u 4 ) + 29Wut(t 2 + u 2 ) + 5444-u 2 £ 2 )M 4 + (17(t 5 + u 5 ) 
-807ut(t 3 + u 3 ) - 2346u 2 t 2 (t + u))M 2 + 2(t + u) 4 (5(t 2 + u 2 ) - ui)] m 10 
+ [ - 768(t + u)M 12 + 800M 14 - 48(i 2 + 32ut + u 2 )M 10 - 48(19(t 3 + u 3 ) - 31tu(< + u))M 8 
+2(787(i 4 + u 4 ) - \22ut(t 2 + u 2 ) - 1546n 2 t 2 )M 6 + (-932(t 5 + u 5 ) + 688ut(t 3 + u 3 ) 
+6004n 2 t 2 (t + u))M 4 + (t + n) 2 (347(t 4 + u 4 ) - 133ut(t 2 + u 2 ) - 592u 2 t 2 )M 2 
-2(< + u) 3 (15(t 4 + u 4 ) + 20ut(t 2 + u 2 ) + 26u 2 t 2 )]m 8 + [256M 16 - 992(i + u)M 14 
+8(209(t 2 + u 2 ) + 34ut)M 12 - 48(55(t 3 + u 3 ) - 29ut(t + u))M 10 + 2(1841(t 4 + u 4 ) 
+704ut(t 2 + u 2 ) - 1722u 2 t 2 )M 8 - 8(395(t 5 + u 5 ) + 478ut(t 3 + u 3 ) - 183u 2 t 2 (t + tt))M 6 
+2(t + u) 2 (782(t 4 + u 4 ) - 808u 2 t 2 - 181ut(t 2 + u 2 ))M 4 - (435(t 7 + u 7 ) + 1327ut(t 5 + u 5 ) 
+1381u 2 t 2 (t 3 + u 3 ) + 617u 3 t 3 (t + u) + 28u 2 t 2 )M 2 + 2(t + u) 4 (21(i 4 + u 4 ) 
+6ui(i 2 + ti 2 ))]m 6 + [32(t 2 - I3ut + u 2 )M 14 - 8(41(i 3 + u 3 ) - 229ut(t + u))M 12 
+8(129(t 4 + u 4 ) - 249ut(t 2 + u 2 ) - QUu 2 t 2 )M 10 - 2(765(t 5 + u 5 ) + 85ut(t 3 + u 3 ) 
-2282u 2 t 2 (t + u))M 8 + 2(t + u) 2 (617(t 4 + u 4 ) - 6Uu 2 t 2 - A8Aut(t 2 + u 2 ))M 6 
-2(285(£ 7 + u 7 ) + 559ut(t 5 + n 5 ) - 315u 2 t 2 (t 3 + u 3 ) - I7f>3u 3 t 3 {t + u))M 4 
+(t + u) 2 (144(t 6 + -u 6 ) + 159ui(i 4 + u 4 ) - 129n 2 t 2 (t 2 + u 2 ) - 280u 3 t 3 )M 2 
-2{t + u) 3 {7{t 6 + u 6 ) + 6ut(t 4 + n 4 ) - u 2 t 2 {t 2 + u 2 ) - 8n 3 i 3 )]m 4 
-tu[- 8(9t 2 - Wut + 9n 2 )M 12 + 8(31(t 3 + u 3 ) + 25ut(t + u))M 10 - 2(141(t 4 + u 4 ) 
+328ut{t 2 + u 2 ) + 334u 2 i 2 )M 8 + 4(25(t 5 + u 5 ) + 84ut(t 3 + u 3 ) + 107n 2 t 2 (i + u))M 6 
+4(t + u) 2 (3(t 4 + u 4 ) + 20ut(t 2 + n 2 ) + 35n 2 t 2 )M 4 - (t + n) 3 (8(t 4 + n 4 ) 
+27ut(t 2 + u 2 ) + 34u 2 t 2 )M 2 + 2(t + u) 4 (£ 2 + ut + n 2 ) 2 ]m 2 
+M 2 tV(i + m) 2 (8M 8 - 16(£ + u)M 6 + 2(t 2 + I2ut + u 2 )M 4 
+ {6{t 3 + u 3 ) - 2ut(t + u))M 2 + tu(t + u) 2 )) 

(5.7) 

For the process g + g ->■ cc( 3 5f ] ) + 66( 3 5[ 8] ), the express of relativistic correction term is too 
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long to give here. We only give the LO cross section. 

. \ / 9 n / „ ,9 . x / „ ,9 \~\ ^ 

X 



F(cc( 3 5l 8] )66( 3 5l 8] )) = a 4 ^ 3 /(108m 3 MV)/ \{ m 2 - t){m 2 - u)(M 2 - t){M 2 - u) 2 



2 



1/ (t + u-2m 2 )(t + u-2M 2 ) x 92M 4 (t + u-2M 2 ) 2 m lb + 4M 2 (767M 8 -U50(t + u)M b 

+18(i9(t 2 + u 2 ) + 99ut)M A -2(86(t 3 + u 3 ) + 279ut(t + u))M 2 + tu(33ut + 7(t 2 + u 2 )))m lA 
+2(2896M 12 - 9812(t + u)M 10 + 4(5417ut + 2899(t 2 + u 2 ) M 8 - 2(2847(t 3 + u 3 ) 
+7877ut{t + u))M 6 + (1003(t 4 + u A ) + 6030u 2 t 2 + 37Wut(u 2 + t 2 ))M A -2tu(19(t 3 + u 3 ) 
+211ut(t + u))M 2 + 76t 2 u 2 (t 2 + ut + u 2 ))m 12 + (3068M 14 - 19624(t + u)M 12 
+(41811(u 2 +t 2 ) + 75706nt)M 10 -9(4421(t 3 +u 3 ) + 11211itt(t + u))M 8 + 2(8729(t 4 + u A ) 
+27487ut{t 2 +u 2 ) + 40854u 2 t 2 )M 6 - (2924(t 5 + u 5 ) + 10065ut(t 3 + u 3 ) + 23123u 2 t 2 (t + u))M A 
+tu(2303ut(t 2 + u 2 ) - 188(t A + u A ) + 3378u 2 t 2 )M 2 - 8t 2 u 2 (84(t 3 + u 3 ) + U9ut(u + t)))m 10 
+(368M 16 - 5800(t + u)M u + 8(2899(t 2 + u 2 ) + 5il7ut)M 12 - 9(4421 (t 3 + u 3 ) 
+H211t**(* + u))Af 10 + (33691(t 4 + u 4 ) + 100610ut(t 2 + u 2 ) + 147662u 2 t 2 )M 8 - 2(6949(t 5 +u 5 ) 
+21495ut(t 3 + u 3 ) + 43712u 2 t 2 (t + u))M 6 + (2274(t 6 + M 6 )+4949ut(t 4 + u 4 ) + 29902u 3 t 3 
+20878n 2 t 2 (t 2 + u 2 ))M 4 + ^(718(t 5 + n 5 )-4033ut(t 3 + u 3 )-8973u 2 t 2 (t + ?i))M 2 
+2i V(679(t 4 + u 4 ) + im9ut{t 2 +u 2 ) + 2164n 2 t 2 ))m 8 - (368(t + u)M lQ - 72(49(t 2 + u 2 ) 
+99ut)M 1A + 4(2847(t 3 + u 3 ) + 7877ut(t + u))M 12 - 2(8729(t 4 + u A ) + 27487ut (t 2 + u 2 ) 
+40854n 2 t 2 ) M 10 + 2(6949(t 5 + u 5 ) + 2U95ut(t 3 + u 3 ) + 43712u 2 t 2 (t + u))M 8 - (5580(t 6 + u 6 ) 
+12481nt(t 4 + u 4 ) + A21Uu 2 t 2 (t 2 + u 2 ) + 58166u 3 t 3 )M 6 + (912(t 7 + u 7 ) - 893ut(t 5 + u 5 ) 
+11979u 2 t 2 (t 3 + u 3 ) + 26506u 3 t 3 (t + u))M A + tu(762(t 6 + u 6 )-5413ut(t A + u A )-21038u 3 t 3 
-15468n 2 t 2 (t 2 + u 2 ))M 2 + 2t 2 u 2 (789(t 5 + u 5 ) + 2413ut(t 3 + n 3 i)3910u 2 t 2 (t + u)))m 6 
+(92(i + w) 2 Af 16 - 8(86(t 3 + u 3 ) + 279ut(t + u)3)M 1A + 2(1003(t 4 + u 4 ) + 6030u 2 i 2 
+3710ut(t 2 + u 2 ))M 12 - (2924(t 5 + u 5 ) + 10065ut(t 3 + u 3 ) + 23123u 2 t 2 (t + u))M w 
+ (2274(t 6 + u 6 ) + 4949ut(t 4 + u 4 ) + 20878u 2 t 2 (t 2 + u 2 ) + 29902u 3 t 3 )M 8 - (912(t 7 + u 7 ) 
-893ut(t 5 + u 5 ) + 11979u 2 t 2 (t 3 + u 3 ) + 26506u 3 t 2 (t + u))M 6 + (152(t 8 + u 8 ) 

- U83ut{t 6 + u 6 ) + 8243u 2 t 2 (t 4 + u 4 ) + 24411n 3 t 3 (t 2 + u 2 ) + 34058n 4 t 4 )M 4 
+tn(334(t 7 + u 7 )-4775ut(t 5 +u 5 )-15750u 2 t 2 (t 3 + u 3 )-26417u 3 t 3 (t + u))M 2 
+2t 2 u 2 (532(t 6 + u 6 ) + 1939ut (t 4 + u 4 ) + 3770n 2 t 2 {t 2 + u 2 ) + 4618u 3 t 3 ))m 4 - tu(-4(7t 2 
+33ut + 7u 2 )M 14 + (76(t 3 + u 3 ) + 8Uut(t + u))M 12 + (188(t 4 + u 4 ) - 2303ut(t 2 + u 2 ) 
-3378u 2 i 2 ) Af 10 + (-718(t 5 + u 5 ) +4033ut(t 3 + u 3 ) + 8973u 2 t 2 (t + u))M 8 + (762(t 6 + u 6 ) 
-5413ut(t A + u A ) - 15468u 2 t 2 (t 2 + u 2 ) - 21038u 3 t 3 ) Af 6 + (-334(t 7 + u 7 ) + 4775ut(t 5 + u 5 ) 
+15750u 2 t 2 (t 3 + M 3 ) + 26417u 3 t 3 (t + n))M 4 + (54(t 8 + u 8 )-2182ut(t 6 + ii 6 )-8161n 2 t 2 (t 4 + ii 4 ) 

- 16101n 3 t 3 (t 2 + u 2 )- !98l2u A t A )M 2 + 54tn(7(t 7 + u 7 ) + 30ut (t 5 + u 5 ) + 68n 2 t 2 (t 3 + n 3 ) 
+99it 3 t 3 (t + n)))m 2 + 2t 2 u 2 (-2M 2 +t + u) 2 (19(t 2 + nt + n 2 )M 8 -(65(t 3 + u 3 ) 

+lllut(t + u))M 6 + (100(t 4 + u 4 ) + 227ut(t 2 + n 2 ) + 300u 2 t 2 )M 4 -27(3(t 5 + u 5 ) + 8ut(t 3 + u 3 ) 
+13u 2 t 2 (t + u))M 2 + 27 {t 2 + ut + u 2 ) 3 )] (5.8) 
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